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Symmetry protected topological (SPT) phases with gapless edge excitations have been shown to 
exist in principle in strongly interacting bosonic/fermionic systems and it is highly desirable to find 
practical systems to realize such phases through numerical calculation. A central question to be ad- 
dressed is how to determine the SPT order in the system given the numerical simulation result while 
no local order parameter can be measured to distinguish the phases from a trivial one. In the tensor 
network approach to simulate strongly interacting systems, the quantum state renormalization algo- 
rithm has been demonstrated to be effective in identifying the intrinsic topological orders. Here we 
show that a modified algorithm can identify SPT orders by extracting the fixed point entanglement 
pattern in the ground state wave function which is essential for the existence of SPT order. The 
key to this approach is to add symmetry protection to the quantum state renormalization process 
and we demonstrate the effectiveness of this algorithm with the example of AKLT states in both 
ID and 2D. 
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Symmetry protected topological (SPT) phases are 
bulk-gapped quantum phases with symmetries pj . If the 
system is on a closed manifold, the ground state does 
not spontaneously break the symmetry. On the other 
hand, if the system has a boundary, there are gapless or 
degenerate edge states as long as the symmetries are not 
explicitly broken. Therefore SPT phases represent a non- 
trivial type of order beyond Landau's symmetry breaking 
theory. Topological insulators and superconductors are 
examples of SPT phases in free fermion systems [2HZ]- 
In one spatial dimension, a complete understanding of 
all possible SPT phases in interacting systems has been 
obtainedJSHHj starting from the classic example of Hal- 
dane phase in spin 1 chains [T5|[T6]. Recently, it has been 
discovered that nontrivial SPT orders can also exist in 
strongly interacting boson/fermion systems in two and 
higher dimensions [TJ [18] . Exactly solvable models were 
presented which has a gapped and symmetric bulk and 
gapless symmetry protected edge states [T5H2TJ] . 

It is highly desirable to find such strongly interacting 
SPT phases in experiments, similar to their free fermion 
counterparts. While the exactly solvable models prove 
the existence of SPT orders in strongly interacting sys- 
tems, it is very unlikely that such models can be realized 
in experiments. In order to determine which physically 
realistic systems can have SPT order, numerical simula- 
tions are necessary. The tensor network renormalization 
algorithm [2TH25] is a powerful and generic approach to 
simulate strongly interacting boson/fermion systems in 
two and higher dimensions and therefore can play a ma- 
jor role in the discovery of strongly interacting SPT or- 
ders. 

A major question to be addressed in the tensor net- 
work approach to simulate SPT phases is how to identify 
the SPT order. Symmetry breaking phases can be iden- 
tified by measuring local order parameters in the ground 
states. However, as SPT ground states do not break any 
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symmetry, no local measurement can distinguish an SPT 
phase from a trivial symmetric phase. Systems with in- 
trinsic topological order (like Z 2 spin liquids or fractional 
quantum Hall systems) on the other hand can be identi- 
fied by measuring the ground state degeneracy on a torus 
[26), [27] or the topological entanglement entropy [281 ES] • 
However, these quantities are both trivial for SPT states. 

An important signature of SPT phases is the existence 
of nontrivial entanglement structure in their ground 
states [U [19]. Compared to trivial symmetric phases 
whose ground states can be simple product states (for 
example the J\(\ t) + I I)) state in the transverse field 
Ising model), the entanglement structure in the ground 
states of SPT phases cannot be totally removed as long as 
symmetry is not broken. Therefore, SPT ground states 
are characterized by short-range entanglement which is 
protected by symmetry. This is similar to systems with 
intrinsic topological orders where the long-range entan- 
glement patterns in the ground states are essential for the 
existence of the order. It has been shown that the long- 
range entanglement patterns can be effectively extracted 
using a quantum state renormalization algorithm based 
on the tensor network representation of the ground states 
[3TH |3"T] . Can similar ideas be applied to SPT orders? 

Naively, one might expect that the quantum state 
renormalizaton algorithm fails to distinguish SPT order 
from trivial symmetric phases, as the algotirhm is de- 
signed to remove short range entanglement structures 
from the state and retain only the long-range one. After 
removing all short range entanglement, the ground state 
of SPT phases becomes a total product state which is the 
same as a trivial symmetry state. However, this is only 
possible if the symmetry of the system is broken during 
the process. If we require that symmetry is always pre- 
served during the renormalization procedure, some short- 
range entanglement structures are always kept which can 
be used to identify the SPT order at the renormalization 
fixed point. 

In this paper, we show how such a symmetry protection 
can be properly added to the quantum state renormal- 
izaiton algorithm. In particular, we demonstrate the ef- 
fectiveness of our algorithm by applying it to the ID and 
2D AKLT phases and show that the SPT order in these 
systems can be successfully identified from the fixed point 
short-range entanglement pattern of the states. The pa- 
per is organized as follows: in section [IIJ we review the 
notion of SPT order (especially that in AKLT state) and 
the quantum state renormalization algorithm which can 
be used to identify intrinsic topological orders; section |H"T| 
discusses how symmetry protection can be added to the 
quantum state renormalization procedure and how it can 
be used to identify the SPT order in ID and 2D AKLT 
states which allows us to determine the phase diagram 
of ID and 2D anti-ferromagnetic spin models; in section 
IV we conclude our discussion and talk about open ques- 
tions. In Appendix [A] we briefly review the notion of 
projective representation; some numerical results of solv- 
ing 2D AKLT-like model are given in Appendix |B| 



II. REVIEW 

A. Symmetry protected topological order 

Symmetry protected topological (SPT) order is char- 
acterized by the robust gapless edge modes of a bulk 
gapped phase without intrinsic topological order or spon- 
taneously symmetry breaking. These gapless edge modes 
are protected as long as the protecting symmetry of the 
system is unbroken. A simple example of an SPT phase 
is the topological insulator in free fermion systems 0- 
4 . For systems with interaction, the SPT phases can 
be classified by group cohomology pQ. The key observa- 
tion used in this classification method is that the edge 
of a d-dimensional SPT phase forms a local d — 1 di- 
mensional system with a special symmetry action. For 
example, in 1-dimensional (ID) SPT phases, the edge 
states are point (O-dimensional) degrees of freedom form- 
ing the projective representations of the symmetry group 
while the physical degrees of freedom in the system all 
form linear representations. Projective representations 
of a symmetry group are classified by the second group 
cohomology H 2 (G,U{\)). Sec Appendix |a| for a brief 
introduction of projective representation. In higher di- 
mensions, symmetry on the edge of an SPT phase can 
act in a non-onsite way related to the higher cohomology 
groups of the symmetry group. Moreover, if translation 
symmetry is present besides on-site symmetry, higher di- 
mensional SPT phases can be obtained by stacking lower 
dimensional ones. In this paper, we are focusing on the 
ID SPT phase with on-site SO(3) rotation symmetry and 
the two dimensional SPT phase with both on-site SO (3) 
symmetry and translation symmetry which can be ob- 
tained by stacking the ID phases along different direc- 
tions. In the following we are going to review our un- 
derstanding of these phases using the example of AKLT 
states. 

One dimensional spin 1 Haldane chain with antiferro- 
magnetic interaction 



H 



Haldane 



= 



(i) 



was the first known example of an SPT phase [T5]. It was 
realized that the bulk of the system is gapped. However if 
the system is put on an open chain, the ground state will 
be four-fold degenerate with two fold degeneracy coming 
from each end of the chain. This result is most easily 
understood at an exactly solvable point in the same phase 
called the AKLT point with Hamiltonian[16 
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The ground state wavefunction was shown to be of a va- 
lence bond structure as shown in FigJTjwhere each spin 1 
decomposes into two spin 1/2 's and the spin 1/2 's form 
singlets |t) — |4)|t) between nearest neighbor sites. The 
total wave function is a product of such singlets between 
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site i and It is then clear that the bulk of the system 
does not break symmetry and is short range entangled. 
Any excitation in the bulk can be gapped. However, if 
the system is open, there will be dangling spin 1/2's left 
at the end of the chain, giving rise to two fold degener- 
acy as long as SO (3) rotation symmetry is not broken. 
Therefore, the AKLT state, as well as the whole Haldane 
phase, has SO (3) symmetry protected topological order. 
It has been shown that the ID Haldane phase is continu- 
ously connected without breaking symmetry to the dimer 
state where each site contains the full Hilbcrt space of two 
spin 1/2's and the spin 1/2's form singlets between near- 
est neighbor sites[32 That is, we can adiabatically 
turn off the projection of two spin 1/2's per site to a total 
spin 1 without affecting the SPT order in the phase. The 
dimer state is a tensor product of local entangled singlet 
pairs and has correlation length. Therefore, it is a fixed 
point representative of the Haldane phase. 




FIG. 1. Valence bond structure of AKLT ground state in 
ID. Each big oval represents a spin 1 and each small circle 
represents a spin 1/2. Connected circles are spin 1/2 singlets. 

The degeneracy of edge spin 1/2 is not only stable 
against any perturbation of the spin 1 chain which is 
SO (3) symmetric, it is also stable even if addition or re- 
moval of local degrees of freedom are allowed as long as 
they form linear representations of SO(3), i.e. they are 
integer spins. It is necessary to consider the second sta- 
bility condition because, when a system is perturbed, it 
is possible that more degrees of freedoms (e.g. orbitals, 
bands) get involved in the interaction. The symmetry re- 
quirement on the system cannot exclude the possibility 
that symmetric local degrees of freedom (linear repre- 
sentations of symmetry groups) to be added or removed. 
The spin 1 Haldane chain is stable even this is allowed as 
the fractional edge spin 1/2 cannot become a trivial spin 
through interaction with any integer spin. However, 
the spin 1 edge state in the spin 2 AKLT state with sim- 
ilar valence bond structure can be removed by adding an 
extra spin 1 to the edge which forms a singlet with the 
original edge spin 1. Therefore, the spin 2 AKLT state is 
in a trivial SPT phase. Its edge state degeneracy is not 
protected under SO (3) rotation symmetry. The differ- 
ence between trivial and nontrivial ID SPT phases can 
be most clearly seen from the edge where the former hosts 
an integer spin and the latter hosts a half integer spin. 
Such a distinction can be generalized to any symmetry 
group where edge states forming projective representa- 
tions correspond to nontrivial SPT phases while edge 
states forming linear representations correspond to trivial 
ones. (See Appendix |A| for definition of projective repre- 
sentation.) The distinction between integer/half integer 
spins (linear and projective representations) is the key to 
our understanding of SO(3) symmetric SPT phases and 



will play an important role in our discussion of quantum 
state renormalizaiton group transforms below. 

The valence bond picture of the ID AKLT state can be 
generalized to 2D [T7] , giving rise to a 2D SPT phase pro- 
tected by both SO (3) and translation symmetry. Con- 
sider the 2D state on a honeycomb lattice with one spin 
3/2 per site. The spin 3/2 can be decomposed into three 
spin 1/2's each forming a spin singlet with another spin 
1/2 on a neighboring site, as shown in Fig. [2] 




FIG. 2. Valence bond structure of AKLT ground state on 2D 
honeycomb lattice. Each big circle represents a spin 3/2 and 
each small circle represents a spin 1/2. Connected circles are 
spin 1/2 singlets. 

The bulk of the system is again symmetric and gapped. 
On the edge of the system, there is a chain of spin 1/2 
left uncoupled. Perturbations to the edge may couple 
these spin 1/2's. However, as long as the perturbation 
preserves SO(3) and certain translation symmetry of the 
honeycomb lattice, the edge state is gapless [10]. There- 
fore, the 2D AKLT state is in a 2D SPT phase protected 
by SO(3) and translation symmetry. The fixed point rep- 
resentative of this phase is again the dimer state which is 
a product of singlet pairs between nearest neighbor sites. 

B. Quantum state renormalization group 

The renormalization group (RG) transformation is a 
mathematical tool to average over the short distance de- 
grees of freedom, and can be used to demonstrate the 
universality of the phases at the RG fixed-point. Extend- 
ing Wilson's original idea of RG, an efficient algorithm 
was proposed, called the density matrix renormalization 
group (DMRG) algorithm [31] to solve for the ground 
state of ID noncritical systems by removing the short 
range entanglement before coarse graining. The DMRG 
algorithm can be conceptually thought of a variational 
algorithm based on matrix product states ansatz. For 
ID critical systems and also systems of higher dimen- 
sion, algorithms based on tensor product states and the 
Multi-scale Entanglement Renormalization Ansatz were 
also proposed to simulate ground state, based on re- 
lated RG ideas of removing local entanglement and coarse 
graining[21 23.. 

Here, we will consider the Quantum State Renormal- 
ization Group (QSRG) transformation which, instead of 
simulating the quantum state of a system, aims to ex- 
tract the universal properties of a system from the de- 
scription of its ground state wave function. The basic 
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idea is to remove nonuniversal short range entanglement 
structures related to the microscopic details of the sys- 
tem from the wave function before coarse graining. After 
many rounds of QSRG the original ground state can flow 
to a simpler fixed-point state. We may then identify or 
even classify the phases using the fixed-point states. Such 
a QSRG algorithm was first shown to be executable based 
on the matrix product state description of ID quantum 
states|SS| which was used for the classification of ID 
gapped phases. The algorithm was further generalized 
to 2D tensor product states, which was demonstrated to 
be effective in identifying symmetry breaking phases and 
phases with intrinsic topological order [3"U]. Our discus- 
sion in this paper is based on these two algorithms and 
we review them below. 

The original procedure of the ID QSRG is as follows 
[3"5] . The algorithm is based on the matrix product state 
(MPS) description of a ID quantum state given by 



i 1 ,i 2 ,....; 



r rr(A il A i *...A i »)\i l i 2 ...i tf ) (3) 



coarse-graining without any truncation for gapped MPS 
with a finite \. 

The matrix U in Q is the aforementioned unitary 
transformation which connects the old and new MPS's 
after one round of QSRG, that is 



A p 



RG 



A' 1 = \ l V l 



(6) 



In summary, after one round of QSRG transformation 
of merging two neighboring sites, the state is trans- 
formed as follows: 



|V) -> |$ = U xa <g> U 3A 



U 2 i-l: 



.#}. (7) 



Then we can combined the new sites and perform another 
round of RG. The procedure is illustrated in FigjS] 
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where ik = l...d with d being the physical dimension of 
a spin at each site, A lk 's are x x X matrices related to 
the physical state \ik) with \ being the inner dimension 
of the MPS. 

In the first step of the QSRG procedure, one merges 
two neighboring sites, says the sites 2i — 1 and 2i into 
a new site. In terms of the matrices A s , that means 
forming a new matrix A as 



a(p<i) — ST^ AP AQ 



(4) 



FIG. 3. Quantum state renormalization group procedure as 
given in [35] . 

For comparison with the 2D version of the QSRG pro- 
cedure, we need to rephrase the steps of the ID QSRG 
procedure using the notion of the double tensor of a MPS. 
The double tensor of an MPS with matrices ^4"s is given 
by 



E a7tf)x = ^Ai x (A 
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One further performs the singular value decomposition 



min(cZ ,x 
1=1 
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where (pq) means the combined physical indices. 

Here, one keeps only the relevant degrees of freedom in- 
dicated by the upper bound min(d 2 , x 2 ) of the above sum. 
These two steps correspond to one round of RG transfor- 
mation. After performing n such steps, The Hilbert space 
of each local site is d 2 -dimensional and corresponds to 
2 n original sites. However, for the noncritical spin chains, 
the block entanglement entropy will be saturated to a fi- 
nite value. This implies that only few degrees of freedom 
couple to the outer sites. From quantum information per- 
spective, there will exist a unitary transformation that 
could transform the d 2 -dimensional space into a smaller 
one. The relevance of the degrees of freedom is indi- 
cated by the magnitude of the eigenvalues of the block 
reduced density matrix. In the MPS representation, the 
dimension of Hilbert space of the block remains bounded 
above by \ 2 as can be seen from the singular value de- 
composition. In this way, it is possible to perform the 



An important property of E is that it uniquely determines 
the matrices, and hence the state, up to a local change 
of basis on each site (HJH2]- That is, if 



E, 



&l,Px 



(9) 



then A l a p and B l a g are related by a unitary transforma- 
tion U: 



B 



a/3 



(10) 



Renormalization transformations on the MPS, including 
the local unitary operations and the coarse graining, can 
be applied through manipulations of E. Treat E as a 
X 2 x x 2 matrix with row index and column index fix- 
Combine the double tensor of the two sites together into 
E = EE. Then think of E a7j( 9 x as a matrix with row 
index a/3 and column index 7%. It is easy to see that 
with such a recombination, E is a positive matrix and 
can be diagonalized 



(11) 
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where we have kept only the non-zero eigenvalues and 
the corresponding eigenvectors Qj g . A is then given by 

= \fm taP , (12) 

which are the matrices representing the new state. Here 
i correspond to the combined index of pq above. Us- 
ing double tensor E, we have implemented the same RG 
procedure for ID MPS as given by Q and (J5J. 

For ID MPS, one can do the exact coarse graining 
without any truncation to the inner dimensions of the 
matrices due to the finite amount of entanglement be- 
tween any segment of the chain and the rest of the sys- 
tem. However, in two dimension we need to find a local 
unitary transformation which removes the short range 
entanglement and discards the irrelevant degrees of free- 
doms before coarse graining. Such a procedure for 2- 
dimensional tensor product state (TPS) has been given 
in |30j . It was shown that the algorithm can remove short 
range entanglement in a 2-dimensional TPS and identify 
the long range entanglement structure, hence the topo- 
logical order, from the fixed point tensors. Here we re- 
view the basic procedure given in |30) . 

A TPS is characterized by an on-site tensor , 
with physical index i and internal indices a/37 etc The 
wave function is given in terms of these tensors by 

|V) = tTr(r i T^..TK..)|M 2 ...z„ l ...) (13) 

i\,ii,...i m ... 

where tTr denotes tensor contraction of all the connected 
inner indices according to the underlying lattice struc- 
ture. Without loss of generality, we consider the honey- 
comb lattice here. The procedure of QSRG is to flow the 
tensor T to its fixed-point form. To do this, we first form 
a positive double tensor Y\ by merging two layers of ten- 
sors T and its hermitian conjugate with the physical 
indices contracted as shown in Fig. [4j namely, 

i 

Note that the resultant double tensor is invariant under 
local unitary transformations by construction. Similarly, 
we can form another double tensor T2 on the neighboring 
site of Ti. Merging Ti and T2 to form a new double 
tensor T with 4 legs, i.e., 

^a'/3'y5',al3-fS = / ] ^^a'-yW ,a-ye^2;e' /3'S' ,e/3S- (15) 
e,e' 

This completes the step (1) of Fig. [4j 

We then spectrally decompose the double tensor, i.e., 
T, i.e., 

T a 'p>YS',alhS = A j(^a'^7'«5')* X "HlW ( 16 ) 
3 

With the above spectral decomposition, we form a new 
tensor T with the physical indices labeled by l,m, n, r, as 




FIG. 4. The QSRG procedure on honeycomb lattice (part 
1): In the step (1), combining double tensor Ti and T2 on 
the neighboring sites into a double tensor T. In the step (2), 
decomposing the double tensor T into the tensor T. In the 
step (3), decomposing the tensor T, resulting in tensors T a , 
T b . 

shown in the step (2) of Fig. |1J i.e., 

j 

The spectral weights encode the relevance of the entan- 
gled components. In this step one has removed the SRE 
by the local unitary transformation which takes the ten- 
sor T into T according to the significance of the spectral 
weights. This is similar to the disentangler in the multi- 
scale entanglement renormalization ansatz (MERA) [37] 
before the coarse-graining. Then, we do a singular value 
decomposition of T in the direction orthogonal to the link 
between Ti and T2 , and decompose tensor T into T a and 
T;, as shown in the step (3) of Fig. |4j i.e., 

T~>lmnr \ rpln ^ rpmr /-t o\ 

"/37<5 — 2-^1 1 a;a,SX X J &;,3eA- \ Lli ) 
A 
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FIG. 5. The QSRG procedure on honeycomb lattice (part 2): 
Merge three tensors T a , Tb and, T c around a triangle to form 
a new tensor T 1 . 

Next, we coarse grain the lattice labelled by the tensor 
T by implementing one step of block decimation in the 
TRG method [21]. To do this, we combine the resultant 
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three tensors that meet at a triangle to form a new tensor 
with physical index / as shown in the Fig. [5j 
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(19) 



where / denotes the combination of all physical indices 
ax, a,2, b\, &2) ci and c%. This new tensor is the resultant 
one after one round of QSRG proposed in [30] . 

After repeating the above QSRG procedure, the origi- 
nal TPS will then flow into a fixed-point state with SRE 
removed but not the LRE. Therefore, we can use this 
method to find out the possible topologically ordered 
phases by examining the fixed-point tensor. Moreover, 
this method is easy to be implemented numerically and 
has been used to identify and study phases with intrinsic 
topological order, for example see [3Tj . 



III. SYMMETRY PROTECTED QUANTUM 
STATE RENORMALIZATION 



Though the above QSRG procedure is quite powerful 
in identifying the intrinsic topological orders, it is not 
suitable when we consider SPT orders. This is because, 
the above QSRG procedure is designed to remove short 
range entanglement, which makes SPT orders indistin- 
guishable from trivial orders. Indeed, consider applying 
the above QSRG procedure to the 2D dimer state on 
the honeycomb lattice. Each bond originally has a dimer 
state on it. When applying the steps in Fig|4j the dimer 
on the horizontal bond shrinks and no dimer state is re- 
generated when splitting in the vertical direction. After 
several rounds of RG, it is easy to see that the dimer 
structure can be completely removed, resulting in a total 
product state with no SPT order. 

FIG. 6. Applying QSRG procedure given in 30 to 2D dimer 
state. After applying steps in Fig(4j the dimer on the middle 
bond gets removed. The dimer structure can be completely 
removed after a couple rounds of RG. 

The key reason that the dimer state flows to a trivial 
state under the QSRG scheme is because the symme- 
try of the system is not preserved when doing local uni- 
tary transformations or coarse graining transformations 
on the state. Therefore, we need to find a way to ex- 
plicitly incorporate the symmetry of the system into the 
RG scheme in order to preserve the short range entangle- 
ment structure of SPT states in the RG flow. In this sec- 
tion, we will propose such a symmetry protected QSRG 
procedure for TPS and demonstrate its effectiveness in 
identifying the SPT order from fixed point tensors. 



In order to devise the 2D procedure, let's start with 
the simpler ID case. The ID QSRG procedure as given 
in [35] does preserve the short range entanglement struc- 
ture in the fixed point tensor of ID AKLT state, hence 
allowing the identification of the SPT order. However, a 
major difference between ID and 2D QSRG is that, in ID 
the amount of entanglement of a segment of an MPS is 
bounded, therefore we can perform coarse graining with- 
out explicitly removing short range entanglement. In 
2D, this is not the case. The amount of entanglement 
of a region in a TPS grows linearly with the length of 
its boundary, therefore straightforward coarse graining 
would lead to unbounded growth of entanglement and 
hence unbounded growth in the size of the representing 
tensor. In order to maintain a bounded numerical com- 
plexity during the RG procedure, it is important that ir- 
relevant entanglement structures are removed using, for 
example, the merging and splitting procedure as shwon in 
Fig. |4j The key to our symmetry protected RG scheme is 
then to preserve symmetry in the merging and splitting 
procedure. Before we do this in 2D, let's first add the 
merging and splitting procedure into ID RG scheme and 
show how symmetry protection can be incorporated. 



A. One-dimensional case 

1. Algorithm 

To construct such a symmetry protected QSRG, the 
key problem is how to preserve the symmetric SRE in the 
process of QSRG. We will adopt the main procedure of 
the ID QSRG as schematically given in Fig. [8] but with 
some essential modifications to preserve the symmetric 
SRE. 

Assume the ground state considered is represented in 
the MPS denoted by a three-rank tensor T % a p, where i is 
the physical index and a , f3 are the left and right indices 
of the inner bond, respectively. In the following we will 
assume translational invariance so that the tensor T is 
the same for all sites. Our proposed symmetry protected 
QSRG is implemented as follows. Firstly, we form a pos- 
itive double tensor by merging two layers of tensors with 
physical indices contracted (as shown in the Fig. [7] after 
step 1). 



(20) 



We then merge two neighboring sites, says the sites i — 1 
and i, into a new block tensor T' (as shown in the Fig. 
[7] after step 3). Perform the spectral decomposition as 
before, i.e., 



(21) 



With such a spectral decomposition we can form a new 
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tensor T as follow: 



T. 



a.b 



E 



(22) 



As mentioned, this step is to perform a local unitary 
transformation on two neighboring sites to remove the 
SRE as the disentangler in MERA. 

Next we need to coarse grain the lattice by first split- 
ting the tensor T 1 into two tensors T-1\ and T- S2 by per- 
forming singular value decomposition. We then merge 
two neighboring T' tensors, say the i-th and the (i + 1)- 
th ones to form a new site tensor T", i.e., 



rplll 

a, 7 



e.7 " 



(23) 



This will complete one round of QSRG as for the 2- 
dimensional case. Schematically the whole procedure on 
the tensors is shown in Fig. [7] and the corresponding 
renormalization operation on the quantum state is shown 
in FigjS] However, the key point is whether we can pre- 
serve the SPT phase during the above process. 
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FIG. 7. 
QSRG. 
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Schematic procedure of the symmetry protected 



To preserve the SPT phase, we should assure that the 
symmetry property for the SPT phase is not destroyed 
during QSRG. Roughly, the above procedure is divided 
into disentangling and coarse-graining. For disentan- 
gling, the on-site symmetry is preserved by the local uni- 
tary transformation as it acts on the double tensor for 
which the physical indices are contracted. 

As for the coarse-graining, it is possible to destroy the 
SPT when splitting and re-merging the neighboring site 
tensors. The symmetry properties of the site tensors for 
a SPT phase will manifest in the degeneracy patterns 
of the singular values when decomposing T into T"s as 



shown in (22) and (23). If there are such degeneracy 



patterns, it means that the SRE encoded in these singular 
values is symmetry protected. We then should be careful 
in assigning appropriate representation for the physical 
legs according to the degeneracy patterns of the singular 
values before re-merging two neighboring sites. To do 
this, we should apply some symmetry generators to the 
MPS to "measure" the quantum number of the site tensor 
and then decide the appropriate representation. 

According to the classification of the SPT phases by 
the using the group cohomology [1] , we should make sure 



^>-<>^-q-^>-<>h>-q- 



^....^....^....^....^ 

| (2) 




J, (3) 

FIG. 8. Schematic procedure of the revised ID QSRG given 



by (20]) -(22]). The step (1) is to remove the SRE by perform- 
ing local unitary transformations. The steps (2) and (3) are 
to coarse grain the lattice by removing the local degrees of 
freedoms and merging two sites. 



that all the site tensors in each step is in the same coho- 
mology class. In practice, this will mean that the physical 
leg of the site tensor should remain in the same cohomol- 
ogy class to preserve the SPT order. Therefore, when we 
split the site tensor T, we should make sure the physical 
leg of the resultant tensor T' is in the same class of the 
group representations as the one for T. If necessary, we 
should enlarge the Hilbert space of the physical leg of T 
when splitting. Moreover, we should assign the appro- 
priate representation for the physical leg by its measured 
quantum numbers so that the tensor T' is in the same 
cohomology class. 

For example, if the on-site symmetry in protecting the 
SPT is 50(3), then the representations for the integer 
spins and the half-integer spins are in different group co- 
homology classes. Therefore, when we split the spin 1 
physical leg of the site tensor T, it may end up with two 
new site tensors whose physical legs can be interpreted 
either as in the spin 1/2 or spin 1 representations. We 
should manipulate the measured J and J z quantum num- 
ber and the possibility of enlarging the Hilbert space to 
make sure that the new site tensors have spin-1 physical 
leg. 

To be more specific about the enlargement of the phys- 
ical Hilbert space, here is the way of doing it: appending 
a singlet state, | Z)fe 1 ,fc 3 =i C kl ,k 2 \k 1 ,k 2 ), where C klM is 
a singlet state coefficient, as schematically shown in Fig. 
Let us denote the state vector for the site tensor T to 
be 

I*) = £ f 0>M> >}|a}|&>| 7 ), (24) 
where a, b together are the physical indices, and a, 7 are 



the inner indices. Then, append a singlet state to enlarge 
the Hilbert space of the physical leg, i.e., 



1*')= E ^,M,7l«)l«)l & )l7) ® ~ E C k u k 2 \kx,k 2 ) 

(25) 



E 9 «,lii,«Afc 2 ,7 l a ) 1°) l fc 2) 1 6) It) • 




M 







(2) 



T* 1 T* H 



a,k b,k 



(6) 
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FIG. 9. A specific way of coarse-graining the MPS in the 
symmetry protected QSRG. 

We can split the new site tensor by performing the 
singular value decomposition so that 

&a,k,a,b,k-f = ^ ^"(o,fc,a),^^^(6,fc,7),/8 (26) 



Then, the new tensors T^t\ and T[ s ^ can be obtained 
from the tensor X and Y by attaching the A, i.e., 



AxI(„, MW , = VA x % fei7)>/3 . (27) 



Finally, we merge two neighboring T' tensors to form a 
new one T" as done in (23). 



2. 1- dimensional AKLT state as the example 

We now demonstrate how to apply the symmetry pro- 
tected QSRG procedure discussed above to a concrete 
example of the SPT phase. Let us consider the 1- 
dimensional AKLT Qi] MPS with T* = a 1 , where a 1 
(i = x,y,z) are the Pauli matrices. This AKLT state 
is shown to have SPT [5J |H1 H3]. Moreover, as shown 
in [35], the fixed-point state of the original QSRG is 
the dimer state characterized by the transfer matrix 
Y,i,j=il\H}(jj\ obtained from its MPS. As we men- 
tioned, this QSRG does not destroy the symmetry and 
so preserves the SPT. Thus, we expect the dimer and 
AKLT states have the same SPT. 

We now apply our symmetry protected QSRG to the 
AKLT MPS. We first perform the disentangler, and then 
the coarse-graining. When coarse-gaining, we find that 
the singular value spectrum of the tensor T has one- and 



three-fold degeneracies. This reflects the AKLT MPS has 
on-site 5*0(3) symmetry protecting its SPT. This degen- 
eracy pattern could be interpreted as combining two spin 
1/2 sites, i.e., 2 (S> 2 = 1 © 3. If we take this interpreta- 
tion and perform the further splitting to obtain two T" 
tensors. Then, the physical leg of T' will be in spin 1/2 
representation. We will obtain the the fixed point double 
tensor, T\'\ = 1, which is a product state. This is in the 
wrong cohomology class from the one of AKLT, whose 
physical leg should be in the integer spin representation. 

To remedy this, we should measure the J and J z 
quantum number of the singular value vectors, and then 
interpret the degeneracy pattern arising from combing 
two spin 1 sites, i.e., 3<8>3 = 1 ® 3 ® 5. We can 
then enlarge the Hilbert space of the physical leg by 
assigning the emptiness to the above 5-representation 
(spin 2) in the enlarged MPS, and identifying the 1- 
representation (spin 0) as l/vo(|l, —1) — 10,0) + | — 1, 1)), 
and the 3-representation (spin 1) as 1/V2(|1,0) — |0, 1)), 
1/V2(|1,-1) - | - 1,1)), and 1/^2(0,-1) - | - 1,0)). 
Here |si,S2)'s with si : 2 = —1,0,1 are the basis vectors 
of 3 Cg) 3. We then identify the labels as the ones for 
the physical legs of the two site tensors T' , respectively. 
In this way, we can split the tensor T to T"'s without 
destroying the SPT and then form the new tensor T" to 
complete one step of our symmetry protected QSRG. 

By this way, we can obtain the fixed point double ten- 
sor T l a J p of AKLT state as follows: 





= 0.524, 


m 


= 0.053, 


^12 


= 0.471, 




= 0.471, 


Til 


= 0.053, 


Til 


= 0.524. 



(28) 

This is the same as the one for the dimer state obtained 
by the same symmetry protected QSRG procedure. It 
then justifies our symmetry protected QSRG in identify- 
ing the SPT. 

We can also perform the coarse-graining of the symme- 
try protected QSRG in the way of Fig. [9] by appending 
a singlet state. In this way, we will obtain the following 
fixed-point double tensor for the AKLT state: 



T}} =0.5, 

Til = OA 



T22 = 0.5, 
T|l = 0-5. 



(29) 

Again, this is the same as the one for the dimer state 
obtained by the same symmetry protected QSRG proce- 
dure. 



Note that the resultant double tensors in ( 28 1 and ( 29 ) 
are different. This reflects the fact that different SREs 
are removed during the process of QSRG. However, once 
we fix the way of QSRG, then both the AKLT and dimer 
states will flow to the same fixed-point double tensor. 



B. Two-dimensional case 

We will now generalize our symmetry protected QSRG 
to the 2D case. We first describe the algorithm for our 
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symmetry protected QSRG on the honeycomb lattice. 
Then, we study a deformed 2D AKLT model on honey- 
comb lattice, and show that there exist a possible SPT 
phase by calculating some order parameters. Finally, we 
demonstrate the power of our symmetry protected QSRG 
by applying it to this AKLT phase. 



1. Algorithm 

The key procedure for the 2D symmetry protected 
QSRG is schematically the same as the one discussed 
in |II B| for non-SPT, and is recapitulated in Fig. [4] and 
Fig. |5| As for the ID case, it contains two steps: dis- 
entangling and coarse-graining. Again, the disentangling 
preserves the on-site symmetries as it acts on the sym- 
metric double tensor. Instead, we should be careful in the 
process of coarse-graining when splitting the tensor T in 
Fig. [4] before re-merging. The symmetry properties of 
the SPT state will manifest in the degeneracy pattern of 
the singular value spectrum when splitting the tensor T. 
We may need to enlarge the Hilbert space for the physi- 
cal legs of the tensor T before splitting, and then split it 
according to the measured quantum number of the phys- 
ical legs. After that, as shown in Fig. [5] we merge the 
three new tensors around a vertex with their physical legs 
in the same cohomology class as the original one. This 
completes one round of the symmetry protected QSRG. 

Again, for 50(3) on-site symmetry the integer spins 
and half-integer spins belong to different cohomology 
classes for classifying SPT phases. Therefore, we should 
assure that the physical legs of all the tensors in the pro- 
cess of QSRG are in the same class of the spin represen- 
tation to keep the SPT. For ID SPT states with 50(3) 
symmetry, physical legs always have to be integer spins. 
However for 2D SPT states with 50(3) and translation 
symmetry, physical legs can be either integer spins or 
half integer spins, depending on the lattice. The key to 
preserving SPT order in the 2D QSRG scheme is to keep 
the spin representation class of the physical leg. Besides, 
the number of degrees of freedom for TPS is much larger 
than the ones for MPS, we should truncate the minor sin- 
gular values to make the numerical computation viable. 
Of course, the truncation should be performed in accor- 
dance with the degeneracy patterns of the singular values 
to ensure the symmetry property of the SPT. Note that 
we are keeping only block translation symmetry (not the 
full translation symmetry) in the RG procedure which is 
sufficient to preserve the SPT order in AKLT states. 



2. 2- dimensional AKLT phase as the example 

To demonstrate the power of our symmetry protected 
QSRG, we should consider a 2D model with SPT phase. 
The simplest one is the 2D AKLT model [T7] on the 
honeycomb lattice, which has 50(3) on-site symmetry. 
Its ground state is the 2D AKLT state which is a pos- 



sible candidate of the SPT phase as expected from its 
ID cousin. Moreover, this AKLT state has a simple TPS 
representation, see [35] for example. 

To be more general, we consider the simple variation 
of the AKLT model with following Hamiltonian 

H = E ^ • % + m£ ■ + • ^) 3 ] ( 3 °) 

<ij> 

This model corresponds to the AKLT one if J\ — 1, J 2 = 
116/243 and J 3 = 16/243. Obviously, this model has 
50(3) on-site symmetry. 



1.0 



AKLT phase 



0.5 



J 0.0- 



-0.5 



-1.0 



-1.0 




- FM phase 



Neel phse 



1.0 



FIG. 10. The phase diagram of the deformed AKLT-like 
model (30 1 by tuning J2 and J3 with Ji = 1 fixed. The labels 
Neel, FM and AKLT stand for the Neel, the ferromagnetic 
and the AKLT phases, respectively. 

We solve this model numerically by using the method 
of simple update [39] for the ground state. We then ap- 
ply the method of TRG [24] to calculate the order pa- 
rameters. A typical phase diagram for the (staggered) 
magnetization is shown in Fig. 
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in which the order 
parameters goes to zero near J3 = 0. This indicates a 
possible SPT phase in this regime. 

After scanning the parameter space of the Hamiltonian 



(30) for the order parameters, we finally obtain its total 
phase diagram as shown in Fig. [lO] It indicates there is a 
phase, we call it the AKLT phase, which does not break 
the on-site and translational symmetry and could be a 
possible SPT phase. Now, we will apply our symmetry 
protected QSRG to give the supporting evidence. 

First, if apply the usual QSRG to the ground state 
of (30), it will flow to the trivial ground state, such as 
This implies that none of the phases in the 
[TO] has intrinsic topological order. 



11 111 — 1 - 

phase diagram Fig. 

Then, we apply the symmetry protected QSRG to the 
2D AKLT state. The TPS of the AKLT state has the 
bond dimension \ — 2. After performing the singular 
value decomposition for the tensor T in Fig. [4J we find 
that only 9 of the 16 singular values are non-zero. More- 
over, these non-zero singular values show the degeneracy 
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FIG. 11. A typical phase diagram for the magnetization M z 
and the stagger magnetization Ml by tuning J2 with J\ = 1 
and J2 = 0.5 fixed. We fix the bond dimension \ — 4 and 
D cut = 24 in this numerical calculation. We also plot the 
difference between the largest two singular values, i.e., Ai — A 2 . 



pattern as 1-, 3- and 5-fold degeneracies. Similar degen- 
eracy pattern appears also for the other ground states in 
the AKLT phase. We also check the ground state for the 
Neel and ferromagnetic phases but find no such degen- 
eracy pattern. This implies that the AKLT phase could 
be a SPT one protected by the on-site SO(3) symmetry 
and the translational invariance. 

Since the original physical leg is in the half-integer 
representation, i.e., spin 3/2, so we should require the 
physical legs of the tensors in the intermediate steps to 
be in the half-integer representation to ensure the SPT. 
The simplest way is to interpret the physical leg of T as 
comprised of four spin 1/2 "partons". Thus, we can un- 
derstand the degeneracy pattern of the singular values as 
follows: 2<g)2®2<g)2 = lffil®3ffl3ffi3ffi5. Accordingly, 
we should enlarge the Hilbert space of the physical leg 
of the tensor T and split it by the rule of the measured 
quantum number in such a decomposition, e.g., fill empty 
elements for two of the three 3-representations and one 
of the two 1-representations. 

By this way of performing symmetry protected QSRG 
on the AKLT state, we finally arrive the following fixed- 
point double tensor: 



T}^ = 0.05; 



npll2 
^112 

n-p221 
Ji 221 

rrp222 

J1 222 



= 0.044; 
= 0.036; 
= 0.05. 



T 



121 = 0.044; 
= 0.036; 



np212 

11 212 



rip 2 11 
^211 

npl22 
1 122 



0.044; 
0.036; 



(31) 



We show that this is also the fixed-point double tensor 
of the 2D dimer state represented by the TPS with the 
tensor T- 



ijk 
ijk 



1 for i,j,k = l, 2. This phase is the valence 
bond solid with each pair of "patrons" on the neighbor- 
ing sites forming singlet state. In the language of TPS, 



the projection operator is the identity. We also apply 
the symmetry protected QSRG to the other numerical 
ground states in the AKLT phases, they all flow to the 
same fixed-point state given by (31). This then demon- 



strate the power of our symmetry protected QSRG and 
also justify the SPT of the AKLT phase. 



IV. CONCLUSION 

After so many years's effort in developing the theo- 
retical perspective of topological order, it is time to ex- 
plore more the possibilities of realizing either intrinsic 
or symmetry protected topological order in real materi- 
als. To motivate the experimental search, it is better for 
theorists to find out the possible realistic models which 
could be fabricated in the real world. As most of the 
exactly solvable higher dimensional models with (sym- 
metry protected) topologically ordered ground states in- 
volve non-realistic interactions, it usually needs numer- 
ical computations to find out the possible topologically 
ordered phases for the more realistic models. 

In this paper we have aimed at such a goal in devising 
a numerical algorithm based on the matrix and tensor 
product state representation of the ground states. Our 
new algorithm of the quantum state renormalization can 
flow a original ground state wave function to its fixed- 
point state in the same class of the SPT. The fixed-point 
state is usually simple and universal due to the removal 
of irrelevant short range entanglements so that it can be 
used to classify the SPT. This algorithm is the modified 
version of the original quantum state renormalization [301 
135] to adopt the symmetry constraints in protecting the 
SPTs. A key observation in our algorithm is to add or 
remove only the symmetry allowed degrees of freedom 
when coarse graining the lattice. 

We have considered the ID and 2D AKLT phases as 
the examples in testing our quantum state renormaliza- 
tion algorithm. These models have on-site SO (3) sym- 
metry so that only integer spins degrees of freedom can 
be added or removed when coarse graining. Our numer- 
ical implantation of the new algorithm on these states 
yield satisfying results by finding that they all flow the 
same fixed-point states as the dimer states do. This con- 
firmed that the AKLT and the dimer states are in the 
same SPT class for the ID and 2D cases. 

Despite of the success of our algorithm, we find that the 
detailed procedure of the symmetry protected quantum 
state renormalization algorithm is not unique, and the 
subtle difference in the procedure will lead to different 
forms of the fixed-point state. Though this is expected 
as the expression of the SRE will depend on the choice 
of the basis and performed local unitary transformations. 
It is then desired to see if there are some criterions to fix 
such ambiguity. 

Finally, it is interesting to apply our algorithm to 
more general realistic models in identifying the new SPT 
phases. Also, it is worthy of devising a new quantum 
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state renormalization algorithm for the phases with both 
the intrinsic and symmetry protected topological orders. 
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parameters and delineate the phase diagrams as shown 
in Fig. 11 and 13 The total phase diagram of this model 



as shown in Fig. TUJ 

We assume the translational invariance for the TPS 
ansatz for the ground state of (30). Using the simple 



update method we can solve the TPS numerically. One 
should caution that the unit cell used in the simple up- 
date is a honeycomb, which is different from the acyclic 
tree of coordination number equal to 3. So, when per- 
forming each step of the simple update, we need to up- 
date the six kinds of tensors and nine bonds as shown in 
Fig. pi 



Appendix A: Projective Representation 

Matrices u{g) form a projective representation of sym- 
metry group G if 

u(gi)u(g 2 ) = w(g 1 ,g 2 )u(g 1 g 2 ), 31,32 £ G. (Al) 

Here ui (51,32) <= U(l) and ui(gi,g 2 ) 7^ 0, which is called 
the factor system of the projective representation. The 
factor system satisfies 



w(ff 2 , ga)u(gi, 32ff3) = u(g 1 ,g 2 )w{gig 2 ,g3) 



(A2) 



for all 31,(72,33 G G. If 0(31,32) = 1, this reduces to the 
usual linear representation of G. 

A different choice of pre-factor for the representation 
matrices it' (3) = (3(g)u(g) will lead to a different factor 
system 0/(31,32): 



^'(31,32) 



ft (3i 32) 
/3(3i)/?(32) 



^(31,32) 



(A3) 



We regard u'(g) and u(g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems 0/(31,32) and 0(31,32) as be- 
longing to the same class ui. 

Suppose that we have one projective representation 
1*1(3) with factor system 01(31,32) of class 01 and an- 
other u 2 (g) with factor system 0*2(31,32) of class 02, ob- 
viously u\(g)®u 2 {g) is a projective presentation with fac- 
tor group ui\ (31, 32)02(31, g 2 ). The corresponding class oj 
can be written as a sum 01 + o 2 . Under such an addi- 
tion rule, the equivalence classes of factor systems form 
an Abelian group, which is called the second cohomology 
group of G and denoted as % 2 [G,U{\)]. The identity 
element 1 £ H 2 [G, U(l)} is the class that corresponds to 
the linear representation of the group. 



Appendix B: Numerical results for solving 
2-dimensional AKLT-like model 



We adopt the method of simple update [39] to solve 
the ground state of ( 30 ) numerically. Then, we apply 
the method of TRG [241 to evaluate the relevant order 




FIG. 12. Diagrammatic representation of the TPS on the 
honeycomb lattice. The tensors with i = 1,2, ..,6 on 

the site labeled by i has three bond indices and one physical 
index. 




FIG. 13. A typical phase diagram for the magnetization M z 
and the stagger magnetization Ml by tuning J3 with Ji = 1 
and J2 = —0.8 fixed. We fix the bond dimension x = 4 and 
D cu t — 24 in this numerical calculation. The left one is with- 
out annealing but the right one is with annealing by tuning 
a small magnetic field in the z-direction when evaluating the 
ground state. This annealing kill the sudden jump of M" . We 
also calculate the energy per site for J3 = —0.2. The ener- 
gies with and without the annealing are —7.303 and —7.150, 
respectively. 

Based on the numerical solution from simple update 
we can further use TRG method to evaluate the expec- 
tation value of the magnetization denoted by M z and 
the staggered magnetization M*. They are the order 
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FIG. 14. 

r = 

are used. 



The spin-spin correlation function as a function of 
. In this numerical calculation, x = 4 and D cll t = 24 



parameters for the ferromagnetic and the Neel phases, 
respectively. In our numerical calculation, we consider 
the bond dimension up to \ — 5 an d keep D cut > x 2 
to ensure the accuracy of the TRG calculation. A typ- 
ical numerical result for M z and M z are shown in Fig. 
[IT] with J\ = 1 and Ji — 0.5, which shows a quantum 
phase transition around J3 = from the Neel phase to 
the AKLT phase as we decrease J3. Note that the dif- 
ference between largest two singular values also indicates 
the quantum phase transition. If we further decrease J3, 
we will reach the ferromagnetic phase. We also plot an- 
other numerical results with J2 = —0.8 in Fig. [13] It 
shows that the Neel order will suddenly drop to zero and 
remain zero after that. On the other hand, the M z is 
zero but suddenly start to grow. However, if anneal the 
result by tuning small magnetic field, then it shows that 
M z grows gradually. This indicates there are only the 
ferromagnetic and the Neel phases in this regime. 

We also calculate the spin-spin correlation function 
the AKLT point in the Fig [l4l It shows the 



at the AKLT point in the Fig [L i 
exponentially-decay behavior as expected for a gapped 
system like the AKLT model. 
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